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Abstract: We study the local linear estimator for the drift coefficient of stochastic differential equations 
^ [ driven by a-stable Levy motions observed at discrete instants letting T — )• oo. Under regular conditions, we 
C/^ [ derive the weak consistency and central limit theorem of the estimator. Compare with Nadaraya- Watson 
i -^H , estimator, the local linear estimator has a bias reduction whether kernel function is symmetric or not under 
j3 ■ different schemes. 
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' 1 Introduction 

= 

■ models in economics and finance, like those for an interest rate or an asset price involve continuous-time 

diff'usion processes. Particularly, their theoretical and empirical applications to finance are quite extensive 
^ ■ (see Jacod and Shiryaev [18]). However, growing evidence shows that stochastic processes with jumps are 

becoming more and more important (see Andersen et al. [2]; Baskshi et al. [7]; Duffie et al. [11]). Recently, 
H ' stochastic processes with jumps as an intension of continuous-path ones have been studied by more and 

more statisticians since the financial phenomena can be better characterized (see Alt-Sahalia and Jacod [1]; 

Bandi and Nguyen [3]). 

A diffusion model with continuous paths is represented by the following stochastic differential equation: 

dXt = ii{Xt)dt + a{Xt)dWu (1.1) 

where Wt is a standard Brownian motion, : R — t- R is an unknown measurable function and a : R — t- R^. 
is an unknown positive function. Many authors have investigated nonparametric estimations for the drift 
function ^(x) and the diffusion function ct(x), which to some extend prevent the misspecification of the model 
(1.1) compare with the parametric estimations. Prakasa Rao [27] constructed a non-parametric estimator 
similar as the Nadaraya- Watson estimator for /u(x). Bandi and Phillips [4] discuss the Nadaraya- Watson 
estimator for these functions of non-stationary recurrent diffusion processes. Fan and Zhang [15] proposed 
local linear estimators for them and obtained bias reduction properties. In a finite sample, Xu [34]extended 
re-weighted idea proposed by Hall and Presnell [16] to estimate <j{x) under recurrence. Xu [33] discussed the 
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empirical likelihood-based inference for nonparametric recurrent diffusions to construct confidence intervals. 
Furthermore, Bandi and Phillips [5] proposed a simple and robust approach to specify a parameter class of 
diffusions and estimate the parameters of interest by minimizing criteria based on the integrated squared 
difference between kernel estimates of the drift and diffusion functions and their parametric counterparts. 

Recently, stochastic processes with jumps have been paid more attention in various applications, for 
instance, financial time series to reflect discontinuity of asset return (see Baskshi et al. [7]; Duffie et al. [11]; 
Johannes [20]; Bandi and Nguyen [5]). In this paper, we consider the stochastic process with jumps through 
the stochastic differential equation driven by an a-stable Levy motion (1 < a < 2): 



where {Zf, t > 0} is a standard a-stable Levy motion defined on a probability space ($7, P) equipped with 
a right continuous and increasing family of cr-algebras {^t,t > 0} and 77 is a random variable independent 
of {Zt}. Z\ has a a- stable distribution S'a(l,^,0) with the characteristic function: 



where fi G [—1, 1] is the skewness parameter. One can refer to Sato [30], Barndorff-Nielsen et al. [6] for more 
detailed properties on stable distributions. Usually, we get observations \Xt-^ti = iA„,i = 0,1,..., n} for 
model (1.2), where A„ is the time frequency for observation and n is the sample size. This paper is devoted 
to the nonparametric estimation of the unknown drift function. Our estimation procedure for model (1.2) 
should be based on {Xt^ , ti = iA„, i = 0, 1, . . . , n}. 

The stochastic differential equation driven by Levy motion has received growing interest from both 
theoreticians and practitioners recently, such as applications to finance, climate dynamics et al.. Masuda 
([23], [24]) proved some probabilistic properties of a multidimensional diffusion processes with jumps and 
provided mild regularity conditions for a multidimensional Ornstein-Uhlenbeck process driven by a general 
Levy process for any initial distribution to be exponential P- mixing. When model (1.2) is specially a 
mean-reverted Ornstein-Uhlenbeck process driven by a Levy process, i.e. fi{x) is known to be linear with 
the form n{x) = 7 — Ax and a = 1, where (7, A) is unknown parameters to be estimated. Based on 
{Xt^,ti = iAn,i = 0, 1, . . . ,n}, Hu and Long [17] studied the least-squares estimator for A > when Z is 
symmetric a-stable and 7 = 0. Masuda [25] considered an approximate self- weighted least absolute deviation 
type estimator for (7, A). Zhou and Yu [37]provcd the asymptotic distributions of the least squares estimator 
of the mean reversion parameter A allowing for nonlincarity in the diffusion function under three sampling 
schemes. However, in model (1.2), the drift function /i(x) is seldom known and the diffusion function may 
be nonlinear in reality. With no prior specified form of the drift function, Long and Qian [22] discussed the 
Nadaraya- Watson estimator for it and obtained the weak consistency and central limit theorem. 

The Nadaraya- Watson estimator given for /x(x) is locally approximating ij,{x) by a constant (a zero- 
degree polynomial). However, in the context of nonparametric estimator with finite-dimensional auxiliary 
variables, local polynomial smoothing has become the golden standard (see Fan [12], Wand and Jones [36]). 
The local polynomial estimator is known to share the simplicity and consistency of the kernel estimators as 
Nadaraya- Watson or Gasser-Miiller estimators but avoids boundary effects, at least when convergence rates 
are concerned. Local polynomial smoothing at a point x fits a polynomial to the pairs {Xi;Yi) for those 
Xi falling in a neighborhood of x determined by a smoothing parameter h. The local polynomial estimator 
has received increasing attention and it has gained acceptance as an attractive method of nonparametric 
estimation function and its derivatives. This smoothing method has become a powerful and useful diagnostic 



dXt = ^i{Xt-)dt + a{Xt-)dZt 



(1.2) 




(1.3) 



tool for data analysis. In particular, the local linear estimator locally fits a polynomial of degree one. In 
this paper, we propose the local linear estimators for drift function in model (1.2). As a nonparametric 
methodology, local polynomial estimator makes use of the observation information to estimate corresponding 
functions not assuming the function form. The estimator is obtained by locally fitting a polynomial of 
degree one to the data via weighted least squares and it shows advantages compared with Nadaraya- Watson 
approach (see Fan and Gijbels [13]). For further motivation and study of the local linear estimator, see Fan 
and Gijbels [14], Ruppert and Wand [29], Stone [32], Cleveland [10]. 

The remainder of this paper is organized as follows. In Section 2, local linear estimator and appropriate 
assumptions for model (1.2) are introduced. In Section 3, we present some technical lemmas and asymptotic 
results. The proofs will be collected in Section 4. 



2 Local Linear Estimator and Assumptions 

We lay out some notations. For simplify, Xi denotes and we shall omit the subscript n in the notation 
if no confusion will be caused. We will use notation "A" to denote "convergence in probability" , notation 
to denote "convergence almost surely" and notation "^" to denote "convergence in distribution" . 

Local polynomial estimator firstly introduced in Fan [12] has been widely used in regression analysis and 
time series analysis. It has gained acceptance as an attractive method of nonparametric estimation of re- 
gression function and its derivatives. The estimator is obtained by locally fitting p-th polynomial to the 
data via weighted least squares and it shows advantages compared with other kernel nonparametric regres- 
sion estimators. The idea of weighted local polynomial regression is the following: under some smoothness 
conditions of the curve m(x), we can expand m{x) in a neighborhood of the point xq as follows: 

m[x) m{XQ) + m [xq)[x - xq) ^ — [x - xq) H 1 j [x - xqY 

V 
j=0 

where = H!^^. 

Thus, the problem of estimating infinite dimensional m(x) is equivalent to estimating the p-dimensional 
parameter Po,Po, - ■ ■ Consider a weighted local polynomial regression: 

n-l p 2 

arg min ^ |y, - ^ ^^-(X, - x)4 K^,^iX,-x), 

where Yi = — - and Kh^{-) = —K{jj-). is kernel function with /i„ the bandwidth. 

What we are interested in is to estimate ft^x) = /3o, hence as Fan and Gijbels [14] remarked, it is 
reasonable for us to discuss p = 1: the local linear estimator for the drift function in this paper. The 
local linear estimator for /Li(x) is the solution Pq of the optimal problem: 

n-l 1 2 



arg min V |l^ - V Pj(,Xi - x)4 Kh^X, 



The solution of /3q is 

^(^) = — — — : ^ ' 

n-1 

where Sn,k = E - x){Xi - x)^ k = 1,2. 

i=0 

We can also write u(x) = 

hn{x) 

Mx)4e^.(x.-.,{||-(^i^)|i}. 
^»(-) - ;i 1 «^(^. -) { ^ - ( V) !f } - 

We now present some assumptions used in this paper. 

(A.l). The drift function /x(-) is twice continuously differentiable with bounded first and second order 
derivatives; the diffusion function cr(-) satisfy a global Lipschitz condition, there exists a positive constant 
C > such that 

\a{y) - a{x)\ < C\y - x\, y, x e M. 

(A. 2). There exist positive constants (Jq and ai such that < uo < cr{x) < a\ for each x G M. 

(A. 3). The solution Xt admits a unique invariant distribution F{x) and is geometrically strong mixing, 
i.e. there exists cq > and p G (0, 1) such that ax{t) < cq/o*, t > 0. 

(A. 4). The density function f{x) of the stationary distribution F{x) is continuously differentiable and 
fix) > 0. 

(A. 5). The kernel function K{-) is nonncgative probability density function with compact support sat- 
isfying: K2 := u^K{u)du < 00, J^^' K'^{u)du < 00. 
(A. 6). As n — > 00, /i — ^ 0, A — > 0, n/iA — > 00. 

Remark2.1. The condition (A.l) ensures that (1.2) admits a unique non-plosive cadlag adapted solution, 
see Jacod and Shiryaev [18]. (A. 3) implies Xt is ergodic and stationary. The mixing property of a stochastic 
process describes the temporal dependence in data. One can refer to Bradley [9] for different kinds of mixing 
properties. For some sufficient conditions which guarantee (A. 3), one can refer to Masuda [24]. The kernel 
function is not necessarily to be symmetric. Sometimes, unilateral kernel function may make predictor easier 
(see Fan and Zhang [15]). 



3 Some Technical Lemmas and Asymptotic Results 

We say that a continuous function G : [0, 00) — >■ [0, 00) grows more slowly than (a > 0) if there exist 
positive constants c, Aoao < « such that G{Xu) < c\"°G{u) for all u > and all A > Aq. 
Lemma 3.1. Let ((){€) be a predictable process satisfying \(j){t)\'^ dt < 00 almost surely for T < 00. We 
assume that cither (p is nonnegative or Z is symmetric. If G{u) grows more slowly than u", then there exist 
positive constants ci and C2 depending only on a, oq, c andXo such that for each T > 

c,E[G{i r m)rdt)^h] < E[Gisup I r 4>is)dZ,\)] < C2^[G(( T |0(t)|"dt)V")]. 

Jo t<T Jo Jo 



Remark 3.1. This lemma can be viewed as a generalization of Theorem 3.2 in Rosonski and Woyczynski 
[28] , where they only dealt with the case that Z is symmetric. 

Lemma 3.2. Suppose that there is a deterministic and nonnegative function $ such that 

l-T 

$"(r) / |(^(t)|"cii 4 1 as T^oo. 
Jo 

Then, we have 

$(T) / |(/)(i)|dZt^5«(l,/3,0). 
Jo 

Remark 3.2. This lemma can be regarded as an extending to a-stable case of Theorem 1.19 in Kutoyants 
[21]. 

Lemma 3.3. Assumptions (A.l) - (A.6) lead to the following result: 

i=0 ^ ^ 

Remark 3.3. In Long and Qian [22], they proved a weaker case: ^ ^h{Xi — x) — > f{x). One can easily 

i=0 

obtain hn{x) — -> K2f'^{x) — {Kif{x))'^ based on this lemma. 

Theorem 3.1. Assume that (A.1)-(A.6) hold and a G (1, 2), then jl{x) ii{x) as n ^ oo. 

Theorem 3.2. Let a G [1,2) and assume that (A.l) - (A.6) are satisfied. 
{i) If {nAh)^~^h = 0(1) and (nA/i)^"« = 0(1) for some k > a, then 

{nAhY-^A{x){i2{x) - n{x)) Sail, 13,0) 

{it) If {nAh)^~»h'^ = 0(1) and (nA/i)^"aA^ = 0(1) for some k > a, then 

(nA/i)i-^ A(x)(A(x) - ii{x) - h^T^ix)) ^ 5«(1, /3, 0) 
where A(x) = \K.-(K,fm^)^-h ^^^^^ ^ (.m^^f^K^K,] _ 

Remark 3.4. In Long and Qian [22], they showed the following results under the Assumptions in this paper: 

(i) If {nAhf'^h = 0(1) and (nA/i)^"« A- = 0(1) for some k > a, then 

(nA/t)^-^ A(x)(A(x) - n{x) - hKi) 5^(1, P, 0) 

(ii) If {nAh)^~ahF' = 0(1) and (raA/i)^~a A^ = 0(1) for some k > a, K{-) is symmetric,, then 

(nA^)i-^ A(x)(A(x) - nix) - h'r^ix)) 5^(1, p, 0) 

where A{x) = — ^, and r^(x) = [/^'(aj)4pr + ^fj," {x)]K2. We can easily observe that the bias 

in the local linear case is smaller than the one in the Nadaraya- Watson case in comparison to the results 
between this paper and Long & Qian [22] whether K{-) is symmetric or not. Furthermore, when a = 1, fi{x) 
is inconsistent with /x(x) easily obtained from Theorem 3.2. 



4 Proofs 



Proof of Lemma 3.1. See Long and Qian (Lemma 2.7). 
Proof of Lemma 3.2. See Long and Qian (Lemma 2.6). 



Proof of Lemma 3.3. 

We first note that 

n-l 



1 " " 
n ^-^ 



X 



1=0 

n-l 



n ^-^ 

i=0 
_ n—\ 

n 



Xi - X 



Xi - X 



fix) / v!'K{u)du 



KhiXi - x) 



^ n—l 

n ^ 



i=0 



Xi - X 



1=0 



Kh{Xi - x) 



Xi - X 



From the stationarity of Xt, we have: 



n— 1 



i=0 



Xi - X 



= E 



f{x) / v!'K{u)du. 



Xi -x ^k 



KhiXi-x)i 



/+0O 
Khiy-x) 
-( 



' — oo 

f+OO 



h 



f{y)dy 



/+0O 
K{u)u''f{x + uh)d 
-oo 

/+00 
u'=i^(u)du. 
-oo 

Thus, from (4.1) and (4.2) it sufHces to prove that 



1 ( 

n-l ( 

Y,\Kh{Xi-x) 

i=0 I 

^ n—l 



n 
1 



Xi - X 



n—l 

-E^ 



i=0 



i^/^(^^ - X) 



Xi - X 



E 



Kh{Xi - x) 



Xi - X 



i=0 



(4.1) 



(4.2) 



(4.3) 



Note that sup |5n,i(a^)| < C^h ^ a.s. for some positive constant Cq < oo by the compact support of 

0<j<n-l 

K{-). Applying Theorem 1.3 (2) in Bosq [8], we have for each integer q G [1, |] and each £ > 



\ n 



n—l 



1=0 



> e 



where 



with p 



4Cn/r 



1\ 1/2 



qaxilp]'^), 



(4.4) 



2q 



and 



+ <^n,[(j+l)p](a;) + ((j + 1)P - [{j + ^)p\)Sn,lij+l)p]+l{x)f- 

Using the Holder inequality and stationarity of Xi , one can easily obtain that s{q) = 0{p'^h~^). 
By choosing q = [\/riKl\/h\ and P= ^ = 0{-\/nh/\fA.), we get 



8^ 



£2 . o{qh) = 0(£^V^). 



(4.5) 



Moreover, we can obtain 



22 1 + 



4Coh- 



1\ 1/2 



qax{\p]A) < C(£)exp(-0(eVnA/i)) 



(4.6) 



under the mixing properties of Xt in (A.3) and (A. 6). 
(4.4), (4.5) and (4.6) imply 



1 



n 



n-l 



i=0 

Therefore, ff: K^iX^-x) (2^Y-lJ^E \Kh{X, - x) (^f 
and (A.6). 



> e < C(£)exp(-0(£2V^^)). 



based on Borel-Cantelli lemma 



Proof of Theorem 3.1. 

It suffices to prove that 

By (1.2), we first note that 

n-l 



i=0 
n— 1 



i=0 

n— 1 



fi=0 



Sn,2 

nh? 

Sn,2 

n. 



h? 



h J nh 
-Xj X \ 



ti + l 



X, ^\Sra^^ I ^^^Xs-)ds + I a{Xs-)dZs) 
U Jt 



h 



^\L{Xi)A+ r^\ti{Xs-) - ii{Xi))ds + a{Xs-)dZ, 



To show the convergence of gn{x), we should prove the following three results: 

(i) gn,iix) [K2f{x) - {Kif{x)f]n{x), as n ^ oo ■ 

(ii) gn,2{x) 0, as n oo ; 

(iii) gn,3{x) 0, as n oo ; 

Proof of (i) : 



1 f 

^ 1=0 ^ 



s, 



n,2 



n-l 



+-S2Kh{Xi-x) 



i=0 



nh? 

Sn,2 



nh^ 



Xi - x \ Snl 

h J nh 



H{x)) 



= : An^i{x) + An^2ix). 

Using Lemma 3.3, it is obviously that 

By the Lipschitz property of fi{x) and the stationarity of Xt, we have 



(4.7) 



(4.8) 



L 



n—l 



'n,2 



n 



h? 



L 



n—l 



+ - V - x\Kh{X, - x) 



n " 

1=0 



Xi-x 



h 



n,l 



nh 



(4.9) 



where L denotes the bound of the first derivative of //(x). 



The two components of the right part are dealt with in the same way, so we only deal with the first one for 
convenience. 



1 n-l ^ ^ n-1 ^ 

- ix, - x\K,{x, - ^) ^ =- E ( - 



1=0 



1=0 

n-l 



n,2 



n 



E 



\Xi - x\Kh{Xi - x) 



n,2 



1 " 
n ^ 



i=0 



Wc find that |Xj - x\Kh{Xi - x) 



Sn,2 



\Xi - x\Kh{Xi - x) 



\Xi - x\Kh{Xi - x) 



S, 



n,2 



nh? 

Sn,2 



(4.10) 



is a.s. uniformly bounded for each i 



by Lemma 3.3 and the compact support of K{-). Similar as the proof of (4.3), we can show that 



n—l 



-y2(\Xi-x\KhiXi-x) 



Sn,: 



n. 



-E 



\Xi - x\Kh{Xi - x) 



Sn,2 



n. 



h? 



0. 



(4.11) 



As for the second part, 



n—l 



hm - V E 



1=0 



\X^ - x\Kh{Xi - X) 
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n 



lim E 



\Xi - x\Kh{Xi - x) 



n,2 



n. 



h? 



= lim K2f{x)E [\Xi - x\Kh{Xi - x)] 



lim hK2f(x ) 



\u\K{u)f{x + uh)du 



/ + 0O 
\u\K{u)du 0. 
-oo 



(4.12) 



It follows that An,2 ^ as n ^ oo. Hence gn,iix) [K2f{x) - {Kif{x)f]n{x) by (4.7)-(4.12). 
Proof of (ii) : 

We first introduce a basic inequality for (1.2): 



/ a{Xs-)dZs 
Jti 

which one can refer to Long & Qian [22], Shimizu & Yoshida [31] and Jacod & Protter [19]. 



sup \Xt -Xt,\< e^^ |/x(Xi)|A + sup 



(4.13) 
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Xi - X 
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Xi - X 



5. 
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nh 

Snl 



\ll{Xs-) - ll{Xi)\d6 

Jti 



nh 
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1=0 
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Sn,2 
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Snl 
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nh 



\Xs--Xi\ds 



sup \Xt-Xt^\ 
ti<t<ti+i 



n 



s 



n,2 



nh"^ 



+ 



Xi - X 



h 



Snl 



nh 



|/x(Xi)|A+ sup 
\ tj<t<tj_i_i 



a{Xs-)dZs 



Similarly as the proof of (4.1), we know that 

_ . ^ 1 



(4.14) 



n 



^KhiXi 
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Sn,2 


+ 


Xi-x 




Snl 
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nK^ 


h 




nh 





\i,{Xi)\-Q{x)AQ, 



where Q{x) = LAe^^\iJL{x)\ 
that is 

n 



K^P{x) + [C^\u\K{u)du)' 

n—l f 

Y,Kn{Xi-x)\ 

i—n ^ 
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Sn,2 
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Xi-x 




Snl 
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nh^ 
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nh 





HXi)\^o. 



(4.15) 



An,4:{x),An,5{x) are dealt with the same approach, hence here we only verify An^/i{x) A 0. 
As for An^4:(x), according to Lemma 3.3, we only need to verify 

n-l 

n 



^ n—l 

-Y,Kh{Xi 



i=0 



x) sup 

ti<t<ti+-i 



I a{Xs-)dZs 
Jti 



0. 



(4.16) 



By Markov inequality and Lemma 3.1, we have 



/ ^ n-l t 

P - V Kh{Xi - x) sup / a{Xs-)dZ^ 
\" ~5 ti<t<ti+i Jti 



n—l 



> e 



< —Y^E sup / Kh{Xi - x)a{Xs-)dZ, 

~0 lh<t<ti+i Jti 

^ / mXi-x)a'^{Xs.)ds) 

n-l 

< ^y2E\Kn{Xi-x)aiA'/- 
ne 'H' L 



i=0 

< O(Ai/")^0. 



(4.17) 



Now, gn,2{x) ^ by (4.14), (4.15) and (4.16). 
Proof of (iii) : 
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By Lemma 3.3, we only need to prove 

n—l 
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nA 
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n-l 



i=0 *» 
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We only prove (4.19) for simplicity. Denote 



Jq" 4>n,l{t,x)dZt . 



i=0 



n—l 



E Kh{Xi-x)Jl'+'a{Xs-)dZ, 

i=0 



we have _ , , „^ , , _ 
With the same argument as the proof of A„_4(a:), we have that: 

— - x) c7{Xs-)dZs >e\=0 ((nA/i)^) . 

We get gnfl{x) ^ by (4.18)-(4.21) and Assumption (A.6). 

Proof of Theorem 3.2. 

Note that 



(nAh)^ aA(x) gn{x) - fj,{x)hn{x) 



{nAhy-^A{x){il{x) - n{x)) 



hn{x) 



hn{x) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



We have obtained hn{x) — > [K2 — Kfjp^x) applying Lemma 3.3, so it is enough to study the asymptotic 
behavior of Bn{x). 

Bn{x) = (nA/i)i-^A(x) gn,i{x) - ii{x)hn{x) + {nAhy-^A{x)gn,2{x) + {nAh)^-^A{x)gn,3{x) 

=: Bn,i{x) + Bn,2{x) + Bn,3{x). (4.22) 

Proof of Bn,i(x) : 

We can express -B„^i(x) as 



n-l ^ 

Bn,i{x) = {nAhf~iK{x)-Y,Kh{Xi - x) I 

^ 1=0 
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Xj — x\ S. 



h J nh 



(M^,)-M^))- 



Using Taylor's expansion, we get 



^i{Xi) - n{x) = ^{x){Xi -x) + -//"(x + OiiXi - x)){Xi - xf, 



where 9i is some random variable satisfying 6i £ [0, 1]. 

Under a simple calculus, we have "^2 ^hi^i - ^) - (^^) th} ' ~ ^) = 0' ^° 



n— 1 
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(4.23) 



By the stationary of Xt, we get 
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One can easily obtain 



Dn,i{x) = i/x"(x)A(x)[(i^2)' - KiKs]{f{x))\nAh)'-^h\l + o(l)). 



Denote 



1 1 """^ 

-Dn,2(x) := -IJ,"{x)A{x)-^^n,i- 



i=0 



(4.25) 



(4.26) 



Note that sup < M\{niS.K) ch a.s. for some positive constant Mi < oo. 

0<i<n-l 

It follows from the proof of (4.3) that - ^ (^n,i — ?• in exponential rate. However, in addition, we should 
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calculate s{q), which may be a little different from s{q). 
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In conclusion, it follows from (4.23)- (4.27) that 
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Similar to the proof of (4.1), under given conditions and Lemma 3.3, we can obtain 
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As in the proof of (4.17) and Lemma 3.3, we can show that 
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It follows from (4.29)- (4.31) and k > a that 

Bn,2ix) = Op{l) ■ {nAh^-^A + Op(l) • {nAh)^-^A^. 

Proof of Bn,3(x) : 
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It follows from the proof of (4.1) that 
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Next we deal with the second term J. By the mean-value theorem, the Lipschitz condition (A.l) and (4.13), 
we have 
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where C is some random variable satisfying G [a{Xs-),a{Xi)] or G [a{Xi),a{Xs-)]. 
Applying Lemma 3.3, and the similar proof of Lemma 3.3, we can obtain 
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hence Ji A 0. As for J2 A- 0, one can refer to (4.31). 
From (4.33)-(4.36) and Lemma 3.2, we have 
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Therefore, it follows from (4.33), (4.37) and Lemma 3.3 that 
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By Slutsky's theorem and remark 3.3, we find 




+ inAh)^~^h'^A{x)T^{x) 



( 




This complete the proof of (ii) in Theorem 3.2. The proof of (i) is similar. 
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